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ABSTRACT 

In this paper we combine archival and proprietary XMM-Newton observations (about 
5deg 2 ) that overlap with the Sloan Digital Sky Survey to explore the nature of the 
moderate- z X-ray population. We focus on X-ray sources with optically red colours 
(g — r > 0.4), which we argue are important for understanding the origin of the X- 
ray background. Firstly, these systems constitute a significant fraction, about 2/3, 
of the z ^ 1 X-ray population to the limit f(2 — 8keV) ~ 2 x 10 _14 ergs _1 cm~ 2 . 
Secondly, their luminosity function under evolution of the form oc (1 + z) 3 suggests 
that they could be responsible for about 17 per cent of the diffuse X-ray background 
to z = 1. Thirdly, their stacked X-ray spectrum in the range 1-8 keV is consistent 
with a power-law distribution with r w 1.4 (without fitting intrinsic absorption), 
i.e. similar to the diffuse X-ray background. We find that the optically red X-ray 
population comprises a mixed bag of objects, both obscured (Nh > 10 22 cm~ 2 ) and 
unobscured (Nh < 10 22 cm~ 2 ), with a wide range of X-ray luminosities up Lx ~ 
10 44 ergs _1 . We argue that dilution of the AGN light by the host galaxy may play a 
role in shaping the continuum optical emission of this population. Finally, we explore 
a possible association of these sources and the moderate- z red (J — Ks > 2 mag) 
AGNs identified in the Two Micron All Sky Survey (2MASS). The median N H of the 
red X-ray sources studied here is « 10 21 cm~ 2 , lower than that found for the 2MASS 
AGNs, suggesting different populations. 
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1 INTRODUCTION 

The origin of the diffuse X-ray background (XRB) remains 
one of the most debated issues of X-ray astronomy. Deep 
X-ray surveys with the new generation X-ray missions, the 
Chandra and the XMM-Newton have revolutionised this field 
demonstrating that the bulk of the XRB can be resolved 
into discrete point sources (e.g. Brandt et al. 2001; Gia- 
conni et al. 2002). These surveys have confirmed previous 
results that luminous (Lx 10 44 ergs _1 ) broad emission- 
line QSOs peaking at redshifts z ~ 1.5 — 2 are undoubtedly 
a major component of the XRB, especially at energies be- 
low « 2keV (e.g. Lehmann et al. 2001). These sources alone 
however, cannot produce the bulk of the XRB. For example 
their steep X-ray spectra (r « 1.9) are inconsistent with the 
spectral shape of the X-ray background (F « 1.4; e.g. Gen- 
dreau et al. 1995) . An additional population of X-ray sources 
is clearly required to account for the XRB properties. 
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There is in particular, accumulating evidence suggest- 
ing that a large fraction of the XRB, particularly at energies 
> 2 keV may arise at redshifts z 1 in moderate luminosity 
systems (Lx Si 10 44 ergs -1 ). For example, previous missions 
at hard energies (> 2keV; HEAOl A-2, Ginga), although 
they lacked imaging capabilities, suggest a statistically sig- 
nificant cross-correlation signal between the detected XRB 
fluctuations and nearby galaxy catalogues (e.g. ESO, UGC, 
IRAS). This implies that about 30 per cent of the XRB could 
be produced at low-z (Jahoda et al. 1991; Lahav et al. 1993; 
Miyaji et al. 1994; Carrera et al. 1995). In the Chandra and 
the XMM-Newton era, a non-negligible fraction of the hard 
X-ray population has been identified with optically extended 
sources suggesting z ^ 1 galaxies (Koekemoer et al. 2002; 
Grogin et al. 2003; Georgantopoulos et al. 2004). This is 
further confirmed by spectroscopic follow-up observations, 
which show a redshift distribution with a peak at z ^ 1 
for hard X-ray selected samples (Barger et al. 2002; Rosati 
et al. 2002; Fiore et al. 2003; Georgantopoulos et al. 2004; 
Georgakakis et al. 2004a). Although the difficulty to spec- 
troscopically identify optically faint sources (R > 25 mag) 
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may skew the above distribution to lower z (e.g. Treister 
et al. 2004), it is nevertheless accepted that these optically 
faint systems (about 25 per cent in the Chandra deep fields) 
cannot drastically modify the observed distribution. 

The evidence above suggests that parallel to the deep 
X-ray surveys targeting the high-z Universe, it is important 
that we also study the nature of the moderate-z X-ray pop- 
ulation. This may indeed hold important clues for the origin 
of the XRB and for interpreting deeper X-ray samples. 

A large fraction of the z 1 hard X-ray sources in par- 
ticular, are associated with optically red galaxies suggesting 
either dust reddening (e.g. Seyfert-2s) or continuum emis- 
sion dominated by stars rather than the central AGN (Geor- 
gantopoulos et al. 2004; Georgakakis et al. 2004a). There 
is also evidence that this population has an average X-ray 
spectrum similar to that of the X-ray background (r w 1.4), 
further underlying its significance for XRB studies (Geor- 
gantopoulos et al. 2004). 

In this paper we explore the nature of the optically red 
nearby (z < 1) X-ray sources and their significance in shap- 
ing the XRB by combining public XMM-Newton data with 
the Sloan Digital Sky Survey (SDSS; York et al. 2000) to 
exploit the uniform 5-band optical photometry and spec- 
troscopy available in this area. Wide field coverage is es- 
sential to probe large enough volume at moderate redshifts 
to provide a representative sample of this class of sources. 
The XMM-Newton with 4 times the Chandra field-of-view 
provides an ideal platform for such a study. Throughout 
this paper we adopt H = 70kms" 1 Mpc" 1 , Qm = 0.3 and 
Q.A = 0.7. 



2 THE X-RAY DATA 

In this paper we use 28 XMM-Newton archival observations 
selected to overlap with the second data release of the SDSS 
(DR2; Stoughton et al. 2002) and with a proprietary pe- 
riod that expired before September 2003. A total of 8 of 
these fields are part of the XMM-Newton/2dF survey (Geor- 
gakakis et al. 2003, 2004b; Georgantopoulos et al. 2004), 
while the remaining 20 pointings are presented by Georgan- 
topoulos & Georgakakis (2005). The exposure times are in 
the range 2-67 ks with a median of about 15 ks. The Galactic 
column density in the direction of these fields varies between 
1.3 - 13 x 10 20 cm" 2 with a median of about 2 x 10 20 cm" 2 . 
A full description of the data reduction, the event file gener- 
ation and the X-ray image production is presented by Geor- 
gantopoulos & Georgakakis (2005). 

For this paper the source extraction is performed in 
the 2-8 keV merged PN+MOS images (when available) us- 
ing the ewavelet task of SAS with a detection threshold of 
6a. This choice of threshold is to minimise spurious detec- 
tions in the final catalogue. The extracted sources for each 
field were visually inspected and spurious detections clearly 
associated with CCD gaps, hot pixels or lying close to the 
edge of the field of view were removed. We further exclude 
from the final catalogue the target of a given XMM-Newton 
pointing (e.g. nearby galaxies or clusters) and a total of 7 
sources that appear extended on the XMM-Newton EPIC 
images and are clearly associated with diffuse cluster emis- 
sion. Fluxes are estimated using an 18 arcsec radius aper- 
ture corresponding to an encircled energy fraction of about 
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Figure 1. Solid angle as a function of 2-8 keV flux (6<r) for our 
survey. 

70 per cent at 1.5 keV. For the spectral energy distribution 
we adopt a power-law with F = 1.8 and Galactic column 
appropriate for each field. Using F = 1.4 instead of F = 1.8, 
typical of the mean spectrum of the XRB and the red AGN 
studied here (see section 4.2), has a minimal effect on the 
estimated fluxes (about 8 per cent) and does not affect the 
results presented here. For the background estimation wc 
use the background maps generated as a by-product of the 
ewavelet task of sas. The final hard X-ray selected sam- 
ple comprises a total of 507 sources above the 6a flux limit 
fx(2 - 8keV) = 5 x 10" 15 ergs" 1 cm" 2 . The area curve 
giving the cumulative area of our survey as a function of 
limiting flux is shown in Figure 1. 



3 THE SAMPLE 
3.1 Selection criteria 

We use the SDSS DR2 catalogue to optically identify the 
hard X-ray selected sources following the method of Downes 
et al. (1986), as described in Georgakakis et al. (2004b) to 
calculate the probability, P, a given candidate is a spurious 
identification. Here we apply an upper limit in the search 
radius, r < 7 arcsec and a cutoff in the probability, P < 0.05, 
to limit the optical identifications to those candidates that 
are least likely to be spurious alignments. 

In this study we use 6 a hard X-ray selected sources with 
2-8 keV flux fx (2- 8keV) > 2 x 10" 14 ergs" 1 cm" 2 . This is 
to ensure sufficient photon-statistics to perform X-ray spec- 
tral analysis for most sources and to minimise the number 
of objects without optical identifications. Out of 230 sources 
above the 6 a detection threshold and fx (2 — 8 keV) > 
2 x 10" 14 ergs" 1 cm" 2 a total of 189 have optical coun- 
terparts to the SDSS limit (82 per cent completeness). 
The choice of the X-ray flux limit above is a trade-off be- 
tween maximum optical identification completeness and suf- 
ficiently large sample size to avoid small number statistics. 
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From the sample above we select sources with red op- 
tical colours, g — r > 0.4 mag, to exclude bluer broad-line 
QSOs. Optically red sources have been identified in non- 
negligible numbers in recent surveys (e.g. Koekemoer et al. 
2002) and may play a key role in shaping the diffuse X-ray 
background. Indeed, this population may be associated with 
obscured AGNs and/or the population of X-ray Bright Opti- 
cally Normal Galaxies (XBONGs; Comastri et al. 2002) that 
are also suggested to harbor deeply buried AGNs (but see 
Georgantopoulos & Georgakakis 2005). Figure 2 plots the 
optical colours of the hard X-ray selected sample with the 
expected tracks for different galaxy types (E/S0, Sbc, Scd, 
Im; Coleman, Wu & Weedman 1980) and optically selected 
QSOs (Cristiani et al. 2004). In this figure it is clear that 
the colour cut g — r > 0.4 mag selects against typical broad 
line QSOs providing a sample that is dominated by sources 
with galaxy-like colours. Nevertheless, high-z QSO (z > 2; 
e.g. Richards et al. 2002; Kitsionas et al. 2005) or reddened 
QSOs (e.g. Wilkes et al. 2002; White et al. 2003) with SEDs 
different to those of optically selected ones may still exist 
within our sample. Such sources are expected to be identi- 
fied by their unresolved (e.g. point-like) optical light profile. 
We also caution the reader that the colour cut g — r — 0.4 
also eliminates from the sample low- (z 0.2) or high-z 
(z J> 0.8) systems with SEDs similar to irregular (Im) type 
galaxies. 

A total of 102 hard X-ray selected sources fulfil the cri- 
teria fx (2-8 keV) > 2xl0~ 14 crgs~ 1 cm~ 2 and 9 -r> 0.4. 
Figure 2 shows that most of them (total of 84) have extended 
optical light profile suggesting moderate-z (^ 1) systems 
where the central AGN does not dominate the optical light 
profile. The 18 sources with unresolved optical light profile 
are associated with either high-z QSOs or Galactic stars. We 
note however, that the SDSS star-galaxy separation is reli- 
able at the 95% confidence limit to r = 21 mag and becomes 
less robust at fainter magnitudes. A total of 25 out of 102 
sources in the sample are fainter that this magnitude limit. 
Also there are only 29 [« 25 per cent, 29/(29+84)] optically 
extended sources with colours bluer than g — r = 0.4 within 
the f x (2 - 8keV) > 2 x 10~ 14 ergs -1 cm~ 2 subsample. The 
red subpopulation, comprising 75 per cent of the optically 
extended X-ray sources, clearly dominates at moderate-z. 

3.2 Redshift estimation 

Spectroscopic redshifts are available for 24 hard X-ray se- 
lected sources with g — r > 0.4 from either the SDSS or our 
own follow-up campaign of the XMM-Newton/2dF survey. 
A total of 22 of the 24 sources are optically extended and are 
assigned moderate redshifts (z < 1). The remaining two X- 
ray sources are optically unresolved and are identified with 
a Galactic star and a z = 1.335 broad-line QSO (Georgan- 
topoulos et al. 2004). For the spectroscopically unidentified 
sources in the sample we use photometric methods to esti- 
mate redshifts. 

In particular for systems with extended optical light 
profile we use the SDSS-DR2 photometric redshifts that are 
based on galaxy templates (Csabai et al. 2002). As a byprod- 
uct of the photometric redshift method each source is as- 
signed a best-fit SED which is a continuous parameter be- 
tween and 1. The two extreme values and 1 correspond 
to ellipticals and actively starforming (irregulars) systems 
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Figure 2. r — i against g — r colour for optically identified 6<r 
sources with fx (2—8 keV) > 2 X 10 -14 erg s _1 cm~ 2 . Filled circles 
and crosses are for optically extended and point-like X-ray sources 
respectively. The expected colours of different galaxy types (E/SO, 
Sbc, Scd, Im) using the template SEDs from Coleman, Wu & 
Weedman (1980) are also shown for redshifts 2 = — 1. Also 
plotted is a QSO template SED for z = — 3 (Cristiani et al. 
2004). Our colour selection g — r > 0.4 is shown with the vertical 
dotted line. 



respectively. Galaxies with intermediate best fit SEDs are 
assigned photometric types between and 1. 

Previous studies on photometric redshifts of X-ray 
sources suggest that galaxy templates work well for sys- 
tems with red optical colours, similar to those studied here 
(Barger et al. 2002, 2003; Mobasher et al. 2004; Georgakakis 
et al. 2004a; Kitsionas et al. 2005). The success rate using 
galaxy templates for these sources is demonstrated in Fig- 
ure 3 which compares the photometric and spectroscopic 
redshift estimates for the optically extended objects with 
available spectroscopic observations. The agreement is good 
with | (zphot — z spec ) |/(1 + z spec ) ~ 0.06. A number of sources 
significantly deviate (<5z > 0.1) from the z spec — z zphot re- 
lation in Figure 3. All of them show broad emission- lines 
and it is likely that AGN light is contributing to their op- 
tical continuum. Figure 3 also shows there is a fair agree- 
ment between the observed spectral classification (emission, 
absorption line) and the galaxy type of the best fit SED 
(early, late types). For absorption- line objects 4 out of 5 are 
best-fit by early-type SEDs (< 0.2). For narrow emission- 
line galaxies and broad- line systems 5/7 and 9/9 are best-fit 
by late type SEDs respectively. One spectroscopically iden- 
tified source has no spectral classification in Figure 3 (see 
below for details). 

A small number of 4 optically faint (r ~ 22 mag; see Ta- 
ble 1) galaxies in our sample are assigned unrealistically low 
photometric redshifts, z fa 0.001, by the SDSS algorithm. 
For these systems we use the relation between X-ray lumi- 



© 0000 RAS, MNRAS 000, 000-000 



4 Georgakakis, Georgantopoulos & Akylas 



0.6 



0.4 



0.2 



I 1 I 1 ^~ 

O Absorption-line spectra^ 
□ Narrow emission-line spectra 
^■Broad-line AGNs 
AUnclassified 
X late-type SED 
• early-type SED 



0.2 



0.4 



0.6 



0.5 



g>0.5 



-1.5 



1 1 


1 i 1 i 1 i \ y 1 
/ / 
/ / 

// 






- 


▲ $5 






I 1 . /l/ 


/ / 
/ 1 

/ i This study 

/ Fiore et al. (2003) 

> 1 > 1 > 1 


40 41 


42 43 44 45 


L x (2 


-8 keV) (erg s" 1 ) 



Figure 3. Photometric against spectroscopic redshift estimates 
for the hard X-ray selected sources with available spectroscopic 
observations. Open circles are systems with absorption line spec- 
tra, open squares correspond to sources with narrow-emission line 
spectra and stars are broad line AGNs. The source in our sample 
(#84 in Table 1) without spectral classification is shown with the 
filled triangle. A cross on top of a symbol indicates a late- type best 
fit SED (photometric type >0.2) to the photometric data, while 
small dots are for early type SEDs (photometric type <0.2).The 
dashed lines are the Sz = ±1 envelope around the z spec = z v ^ ot 
continuous line. 



Figure 4. 2-8 keV X-ray luminosity against X-ray— to-optical flux 
ratio for red X-ray sources with available optical spectroscopy. 
Open circles are systems with absorption line spectra, open 
squares correspond to sources with narrow-emission line spectra 
and stars are broad line AGNs. The source in our sample (#84 
in Table 1) without spectral classification is shown with the filled 
triangle. The continuous line is the best-fit relation to the data. 
The dashed line is the best-fit relation of Fiore et al. (2003) for 
their sample of non type-f AGNs. 



nosity and X-ray-to-optical flux ratio (see next section for 
the definition of log fx / '/opt) suggested by Fiore et al. (2003) 
to get a rough redshift estimate. Figure 4 plots this relation 
for the optically red sources in our sample with spectroscopic 
redshift available. Using this relation we estimate z « 1 — 2 
for the 4 optically faint sources. These systems are marked 
in Table 1. The redshift uncertainty using the la rms scatter 
around the best-fit relation in Figure 4 is estimated to be 
Sz/(l + z) «0.2. 

Figure 5 plots the redshift and spectral type distribu- 
tion of the optically extended subsample. About half of the 
sources lie at z <, 0.4 and have best fit spectral types < 0.4. 

A total of 18 X-ray sources with red colours have point- 
like optical light profile and are most likely associated with 
either Galactic stars, or high- z QSOs. Using the photometric 
methods described by Kitsionas et al. (2005) we find that 
5 of these sources are best fit by stellar templates with the 
remaining assigned photometric redshifts z > 1.3. These 18 
sources will not be considered in the rest of the paper. 

We also exclude from the analysis source #47 in Table 
1 at 2 = 0.005 with L x ~ 2 x 10 39 ergs _1 . This system 
has X-ray-to-optical flux ratio < — 2 (sec Figure 6 below) 
and has been classified 'normal' galaxy by Georgantopou- 
los et al. (2005) with X-ray emission associated with stellar 
processes. For completeness this is included in Table 1 (see 
below). The final sample used in the analysis comprises a 
total of 83 systems with g — r > 0.4 and extended optical 
light profile. For luminosity estimates photometric redshifts 
are used unless spectroscopic redshifts are available 
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Figure 5. Redshift (photometric or spectroscopic) and photomet- 
ric type distribution for the red optically extended X-ray sources. 

3.3 X-ray spectral analysis 

We explore the X-ray spectral properties of our sample us- 
ing the XSPEC vll.2 package. For sources with small num- 
ber of net counts we use the C-statistic technique (Cash 
1979) specifically developed to extract information from low 
signal-to-noise ratio spectra. The data are grouped to have 
at least one count per bin. We adopt an absorbed power- 
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law model (wabs*pow) to constrain the absorbing column 
density Nh by fixing the power law index to F — 1.8, i.e. 
in-between radio loud and radio quiet AGNs. For sources 
with sufficient net counts (about 150) we perform standard 
X 2 spectral fitting. The data are grouped to have a mini- 
mum of 15 counts per bin to ensure that Gaussian statistics 
apply and we require that the spectrum has at least 10 spec- 
tral bins. In the case of \ 2 analysis an absorbed power-law 
(wabs*pow) is fit to the data yielding both the N H and 
the power-law photon index V. For both the \ 2 an d the 
C-statistic analysis the fit was performed in the 0.2-8 keV 
energy range where the sensitivity of the XMM-Newton is 
the highest. The estimated errors correspond to the 90 per 
cent confidence level. 

The column densities estimated above are in the ob- 
server's frame and therefore lower than the rest-frame Nh 
because the fc-effect shifts the absorption turnover to lower 
energies. The relation between intrinsic and observed col- 
umn density scales approximately as (1 + z) 2,65 (e.g. Barger 
ct al. 2002). This redshift correction is applied to all sources, 
after subtracting from the observed Nh the appropriate 
Galactic column in the direction of the source. The rest- 
frame Nh is then used to estimate the unabsorbed X-ray 
luminosities. 

The optical and X-ray properties of the hard X-ray 
selected sample with g — r > 0.4 and /x(2 — 8keV) > 
2 x 10 -14 ergs -1 cm -2 used in this paper is presented in 
Table 1 which has the following format: 

I. Identification number. 

2-3. Right ascension and declination of the X-ray cen- 
troid position in J2000. 

4. SDSS r-band magnitude of the optical counterpart. 

5. Probability, P, the optical counterpart is a chance 
coincidence. 

6. 2-8 keV X-ray flux in ergs -1 cm -2 . 

7. X-ray-to-optical flux ratio, log/x// op t, estimated 
from the relation 

log ^- = log f x (2 - 8 keV) + 0.4 r + 5.60. (1) 

J opt 

The equation above is derived from the X-ray-to-optical 
flux ratio definition of Stocke et al. (1991) that involved 0.3- 
3.5 keV flux and V-band magnitude. These quantities are 
converted to 2-8 keV flux and r-band magnitude adopting a 
mean colour V — R = 0.7, the transformations of Fukugita 
et al. (1996) and a power-law X-ray spectral energy distri- 
bution with index T = 1.8. 

8. Spectroscopic redshift measurement if available. 

9. Spectroscopic classification. AB is for absorption line 
systems, NL is for narrow emission-line sources and BL sig- 
nifies broad emission-line optical spectra. Source #84 has a 
spectroscopic redshift estimate from the NED but no spec- 
tral classification is available. 

10. Photometric redshift estimate. 

II. 2-8 keV X-ray luminosity in units of ergs -1 , cor- 
rected for X-ray absorption. 

12. Rest-frame column density, Nh, estimated from the 
spectral fitting analysis in units of 10 22 cm -2 . 

13. Power-law spectral index, T, for those sources with 
sufficient net counts (about 150) to perform \ 2 analysis. The 
median T for these sources is « 1.8. 
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Figure 6. r-band magnitude against 2-8 keV flux for the opti- 
cally extended g — r > 0.4 sources. The lines log fx / j opt = ±1 
delineate the region of the parameter space occupied by powerful 
unobscurcd AGNs. 
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Figure 7. Nh versus photometric redshift spectral type. Ellip- 
tical SEDs correspond to spectral type=0 while irregulars are 
assigned type=l. Upper limits in the hydrogen column density 
Nh are shown with an arrow. A cross on top of a symbol is for 
sources with broad-line optical spectra. Similarly, open squares 
and open circles on top of a dot correspond to sources with nar- 
row emission-line and absorption optical spectra respectively. 
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Table 1. Optical and X-ray properties of the hard X-ray selected sample. 



4 RESULTS 

4.1 Optical and X-ray properties 

Figure 6 plots r-band magnitude against 2-8 keV X-ray flux. 
The log(/x//o P t) = ±1 lines in this figure delineate the 
region of the parameter space occupied by AGNs (Stocke et 
al. 1991). In this figure the red hard X-ray selected sources 
have a wide range of X-ray-to-optical flux ratios. This may 
suggest diverse AGN accretion rates, black hole masses, line- 
of-sight obscurations and host galaxy properties. 

The present hard X-ray selected sample is indeed, het- 
erogeneous in its optical and X-ray properties. Figure 7 plots 
rest-frame Nh as a function of best-fit SED estimated as a 
by-product of the photometric redshift estimation. The X- 
ray sources studied in this paper are hosted by a wide range 
of galaxy types with no strong correlation with intrinsic X- 
ray absorption. Figure 8 demonstrates that the unobscured 
2-8 keV X-ray luminosities, providing a proxy to accretion 



rate and black hole mass, take a wide range of values and 
are uncorrelated with the estimated Nh- 

We further explore the nature of the red optically ex- 
tended hard X-ray selected sources using the optical spec- 
troscopy available for a small but representative subsam- 
ple (21 out of 83 sources). Most of the narrow emission- 
line galaxies are associated with obscured systems (Nh 
10 21 cm -2 ; see Figures 7, 8), suggesting Seyfert-2 type activ- 
ity. It is likely that the optical emission of these systems is 
dominated by the host galaxy rather than the central engine 
(Barger et et al. 2002; Georgakakis et al. 2004a). One of the 
narrow emission-line galaxies nevertheless, does not show ev- 
idence for X-ray obscuration (Nh ?S 5 x 10 20 cm~ 2 ) although 
the X-ray luminosity of L x (2 - 8keV) « 2 x lO^ergs" 1 
is higher than the maximum that starburst activity is be- 
lieved to be able to produce (« 10 42 ergs _1 ; Zezas, Geor- 
gantopoulos & Ward 1998; Moran et al. 1999). A similar 
class of sources are the 'composites' discussed by Moran et 
al. (2002) and Georgantopoulos, Zezas & Ward (2003) where 
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photometric redshift from the Lx — log fx / fopt relation 
2 source #72 is classified Galactic star by 2QZ. Inspection of the optical spectrum however, suggests an early type galaxy at z — 0.243 
3 sourcc #84 has no spectral classification in NED 



star-formation is suggested to outshine the AGN. Although 
this is a possible scenario we cannot use diagnostic emission 
line ratios for classification because the H/3 line is diluted by 
instrumental features. 

Sources with absorption-line optical spectra have 
Lx(2 — 8keV) > 5 x 10 42 ergs _1 and are likely to be- 
long to the class of X-ray Bright Optically Normal Galax- 
ies (XBONGs; Comastri et al. 2002; Georgantopoulos & 
Georgakakis 2005). None of these sources show evidence for 
large obscuring columns. This is in agreement with Geor- 
gantopoulos & Georgakakis (2005) who claim that this pop- 
ulation comprises a significant fraction of sources where the 
AGN emission is diluted by the host galaxy. We neverthe- 
less find a number of X-ray sources with early-type best- 
fit SED (type< 0.2) and high obscuring column densities 
(Nh > 10 22 cm~ 2 ; see Figure 7). Although optical spec- 
troscopy is not available for these systems they are candi- 



dates for the class of obscured XBONGs similar to the 'Fiore 
P3' source (Fiore et al. 2000). 

Within the spectroscopic subsample there are a num- 
ber of low- z sources with unabsorbed X-ray spectra, broad 
emission-lines and optical colours somewhat redder than 
those of typical QSOs. This suggests either dust reddening 
and/or optical continuum dominated by the host galaxy. For 
7 of these sources flux calibrated optical spectra are available 
allowing estimates of the optical reddening of the broad-line 
region using the Ha/H/3 emission line ratio. Table 2 lists 
the resulting E(B — V) along with the corresponding Nh 
assuming a Galactic dust-to-gas ratio (Bohlin, Savage & 
Drake 1978). For some sources the column density listed in 
this table is higher than the upper limit estimated from the 
X-ray spectrum. This may suggest either anomalous dust- 
to-gas ratios (Maiolino et al. 2001; Maiolino, Marconi & 
Oliva 2001; Georgantopoulos et al. 2003) or the presence of 
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Figure 8. Nu versus unabsorbcd 2-8 keV X-ray luminosity. The 
symbols are the same as in Figure 7. 



a warm-absorber (Crenshaw, Kraemer & George 2003). We 
note however, that stellar absorption from the host galaxy 
may affect the measured Balmer emission line fluxes, par- 
ticularly the HP, overestimating the optical reddening and 
hence, the corresponding Nh- Nevertheless, only 2 of the 7 
broad-line sources show evidence for large reddening. This 
suggests that dust extinction alone cannot fully explain the 
red colours of all the spectroscopically identified broad line 
systems. 

We explore the possibility of stellar light contamination 
of the optical colours using the prescription below: for the 
intrinsic spectrum of the central engine we adopt the SDSS 
composite QSO SED (Vanden Berk et al. 2001). We then 
estimate the contribution of the AGN to the broad band 
magnitudes under the conservative assumption that the ob- 
served u-band flux is entirely due to AGN emission. We shift 
the SDSS composite QSO SED by the appropriate redshift 
and redden it by the observed extinction listed in Table 2, 
adopting the Galactic extinction law. This spectrum is then 
scaled so that the integrated u-band flux matches the ob- 
served one. Finally, for the scaled spectrum we estimate the 
flux through the SDSS r-band and compare it with the ob- 
served one. The estimated fraction of QSO emission in the 
r-band as well as the corresponding aox^ indices are listed 
in Table 2. We note that both the resulting flux ratios and 
the aox are conservative upper limits since we assume that 
all the observed it-band flux originates from the central en- 
gine. With the exception of the two most obscured sources, 
the AGN upper-limit contribution to the r-band is in the 
range 40 — 60 per cent. For the two sources with the highest 
E(B — V) estimate in Table 2 the estimated aox upper lim- 

t we define a Q x = log(/250o//2k c v)/2.605, where /2500 and 
fl kcV is the fl ux density at 2500A and 2 keV respectively. 



Figure 9. g — r colour as a function of redshift for spectroscopi- 
cally identified broad-line systems in our sample. The curves are 
the colour- z tracks of galaxy/QSO composite spectra with frac- 
tional contribution of the galaxy emission at 6000A fixed to 80 
per cent. Different galaxy types are plotted: the dotted line is for 
E/S0, the short-dashed line corresponds to Sbc, the long dashed 
and dot-dashed lines are for Scd and Im types respectively. The 
upper and lower continuous lines represent pure E/S0 and QSO 
spectra and are plotted for comparison. 



its are unrealistically high and do not place any constraints. 
Adopting a more realistic aox index of 1.5 (Brandt, Laor 
& Wills 2000) for these two systems we estimate an AGN 
contribution of S3 1 per cent in the r-band. This increases 
to an upper limit of ft* 10 per cent in the extreme case of 
aox = 2.2, i.e. typical of obscured systems (e.g. BAL QSOs; 
Brandt et al. 2000). 

Figure 9 shows how the g — r colours of the SDSS QSO 
spectrum are modified at different redshifts by superimpos- 
ing galaxy SEDs spanning the range E/SO-Im at a fixed 
galaxy/AGN r-band light ratio of 80 per cent. Templates 
with a significant old stellar component (e.g. E/S0, Sbc) can 
reproduce the observed g — r colours of most sources. Late- 
type galaxy templates (e.g. Scd, Im) produce composite 
spectra that are consistent with the colours of some sources 
in Figure 9. These late-type SEDs however, require on av- 
erage larger galaxy/AGN light fractions (> 80 per cent) to 
reproduce the mean g — r colours of the sources in Figure 
9. We conclude that dilution of the optical AGN emission 
by stellar light is a possible explanation for the observed red 
colours of moderate- z broad-line systems if the galaxy/AGN 
ratio is larger than « 80 per cent in the r-band. 



4.2 Mean X-ray spectral properties 

We attempt to constrain the co-added spectral properties of 
the red optically extended sources to explore whether they 
are consistent with those of the XRB (r w 1.4; Lumb et 
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ID 


Ha/H/9 


E(B - V) 


Nx 1 
(cm- 2 ) 




r-band 
AGN fraction 




o n 

z.y 


U.UU 




<1. i z 


<U.04 


66 


3.6 


0.13 


6.5 x 10 2() 


<1.52 


<0.46 


06 


2.7 


0.00 




<1.14 


<0.40 


08 


4.7 


0.34 


1.7 x 10 21 


<1.67 


<0.56 


45 


18.2 


1.47 


7.4 x 10 21 


<2.77 


<3.36* 


51 


20.9 


1.58 


7.9 x 10 21 


<2.98 


<12.01* 


83 


2.8 


0.00 




<1.40 


<0.60 



The Njj listed here is estimated from E(B — V) adopting the Galactic 
dust— to-gas ratio (Bohlin, Savage & Drake 1978) 

*For these two systems the AGN emission in the r-band is estimated, 
to be higher than the observed flux. 



Table 2. Optical and X-ray properties of the 7 broad emission-line sources with flux calibrated spectra 



XMM-Newton 


number of 


r 


xl 


d.o.f. 


filter 


sources 








THIN 


56 


1 4 7 +0.04 
i - 4 '-0.04 


1.08 


562 


MEDIUM 


27 


1 40+ - 05 


1.06 


335 



Table 3. The mean X-ray spectral properties of the red optically 
extended hard X-ray selected sample. Sources observed through 
the THIN and MEDIUM XMM-Newton filters are merged sepa- 
rately. The errors correspond to 90 per cent confidence level. 



al. 2002; Kushino et al. 2002). The individual spectra of 
the sample are merged using the mathpha task of FTOOLS 
to produce 3 independent coadded spectral files for the 
PN, MOS1 and MOS2 detectors respectively. Sources ob- 
served through the MEDIUM and THIN filters of the XMM- 
Newton are merged separately. The combined spectra arc 
grouped to a minimum of 15 counts per bin to ensure that 
Gaussian statistics apply. The auxiliary files of individual 
sources were combined using the addarf task of ftools. 
Using the xspec vll.2 software, we fit a single power-law to 
the data absorbed by the Galactic column of 2 x 10 20 cm~ 2 
(wabs*pow) . The fit is performed in the 1-8 keV range. This 
is to avoid energies below 1 keV where the XRB spectrum is 
likely to change shape (McCammon & Sanders 1990; Gen- 
dreau et al. 1995). The results for both the THIN (see also 
Figure 10) and the MEDIUM filters respectively are pre- 
sented in Table 3. The best-fit power-law index of the two 
separate set of spectra are in fair agreement within the 90 
per cent confidence level. The mean X-ray spectrum of the 
population is consistent with Y ~ 1.4 similar to that of the 
diffuse XRB underlying the importance of this population. 



4.3 X-ray luminosity function 

In this section we derive the binned X-ray luminosity func- 
tion of the z < 0.4 subsample using the standard non- 
parametric \ IV max method (Schmidt 1968). The choice of 
redshift cutoff is firstly to explore the low-z luminosity func- 
tion of these systems and secondly to minimise the incom- 
pleteness of our sample due to optically unidentified sources, 
many of which are expected to lie at z J; 1. For a given red- 





1 2 5 

channel energy (keV) 

Figure 10. The EPIC-PN mean X-ray spectrum of the optically 
red hard X-ray selected sample. Only sources observed through 
the THIN XMM-Newton filter are presented here. The continuous 
line is the single absorbed power-law (wabs*pow) model fit to the 
data. 



shift and X-ray luminosity interval the binned luminosity 
function is estimated from the relation: 



>dlogL = ^ 



Vrn 



i{L,z)' 



(2) 



where V ma x.i(L, z) is the maximum survey volume available 
to the source i of the sample. This is estimated from the 
relation 



C 

~Wo 



r 

J zl 



n(L, 



dV 

dz 



dz dL, 



(3) 



where dV/dz is the volume element per redshift interval dz. 
The integration limits zl, z2 are the minimum and max- 
imum redshifts possible for a source of luminosity Lx to 
remain within the flux limits of the survey and to lie within 
the redshift bin. Q(L, z) is the solid angle of the X-ray sur- 
vey available to a source with luminosity L at a redshift 
z (corresponding to a given flux in the X-ray area curve). 
The logarithmic bin size of the luminosity function varies 
so that each bin comprises approximately equal number of 
sources, N « 10. The uncertainty of a given luminosity bin 
is estimated from the relation: 
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E 



Vmax,i Z) 



(4) 



In addition to the non-parametric method above we also de- 
rive the luminosity function using the parametric maximum 
likelihood method (Tammann, Sandage & Yahil 1979) as- 
suming a conventional double power law (e.g. Barger et al. 
2005; Sazonov & Revnivtsev 2004) of the form: 



d<f>(Lx, z) _ ^* 

dlogLx ~ (£/£*) 91 + (L/L+y 2 ' 



(5) 



where the characteristic luminosity evolves with redshift ac- 
cording to the relation = L (1 + z) A . Our sample only 
covers a narrow redshift range that does not allow constrain- 
ing of the evolution parameter A. We therefore, choose to 
fix A — 3, i.e. similar to that estimated by Barger et al. 
(2005) and Ueda et al. (2003). We also assume that 0*, the 
normalisation of the luminosity function, is constant with 
redshift (e.g. no density evolution; Barger et al. 2005). This 
is then estimated from the relation 



<t>* = Ngal/ 



fi(L, z) dL dV/dz dz 



(6) 



where N gai is the total number of galaxies in the survey and 
fl(.L, z) is the solid angle of the X-ray survey available to a 
source with luminosity L at a redshift z, i.e. the area curve 
at different flux limits. The best-fit parameters are shown 
in Table 4. The uncertainties correspond to the 68 per cent 
confidence level and are estimated by fixing one parameter 
at a time and then finding the maximum likelihood fit for 
the remaining two. This is repeated by varying the value of 
the fixed parameter between successive runs. The la errors 
for that parameter correspond to SL — 2.7/2 regions around 
the maximum likelihood fit. We note that the 68 per cent 
upper limit of the power-law index is unconstrained. This 
is because our sample does not comprise large number of 
luminous systems (Lx > 10 44 ergs _1 ) to provide tight upper 
limits on <?2- Larger area coverage is required to provide 
large numbers of such systems at z < 0.4. We also note that 
heavily obscured AGNs (Nh ii 10 23 cm~ 2 ) are likely to be 
underrepresented in our 2-8 keV selected sample. 

The results of our analysis are shown in Figure 1 1 along 
with the maximum likelihood fit to the data. Also shown 
in this plot are the X-ray luminosity functions of Barger 
et al. (2005) for their 0.1 - 0.4 redshift bin and the local 
X-ray luminosity function of Sazonov & Revnivtsev (2004) 
shifted to the 2-8 keV band assuming T = 1.8 and evolved 
to z = 0.235, the median redshift of our subsample, assum- 
ing luminosity evolution of the form oc (1 + z) 3 . There is 
fair agreement between our luminosity function and those 
derived in previous studies. This suggests that the red sub- 
sample comprises a sizable fraction of the X-ray population 
at z < 0.4. In section 3.1 we indeed find that about 75 per 
cent of the X-ray sources with optically extended counter- 
parts (most likely to lie at low-z) have red optical colours. 
This does not suggest that broad-line AGNs, which typically 
have blue optical colours, are an unimportant component of 
the X-ray population. Indeed, the colour selected subsample 
studied here also comprises unboscured broad-line AGNs, 
which appear optically redder most likely because of stellar 
light contributing to the broad-band colours. The evidence 
above underlines the importance of the red sub-population 




10 42 1Q 43 

L x (2-8 keV) 



10 44 
(erg s" 1 ) 



Figure 11. The luminosity function (filled circles plotted at the 
median Lx of each bin) for the optically red X-ray population 
with redshifts z < 0.4. The continuous line is the maximum like- 
lihood fit to the data. This is compared to the Barger et al. (2005) 
X-ray luminosity function for all their sources in the range 0.1—0.4 
(open circles). The dashed line is the maximum likelihood fit to 
the Barger et al. (2005) data. The dotted line is the local RXTE 
luminosity function estimated by Sazonov & Revnivtsev (2004). 
This is shifted to the 2-8 keV band assuming T = 1.8 and evolved 
to z = 0.235, the median redshift of our subsample, assuming 
luminosity evolution of the form oc (1 + z) 3 . 



for understanding the different types of sources that make 
up the diffuse X-ray background. 

Using the best-fit parametric form for the luminos- 
ity function we derive the X-ray emissivity (luminosity per 
Mpc 3 ) 



L dL, 



(7) 



also taking into account the column density distribution of 
our sources. We also estimate the fractional contribution to 
the 2-8 keV X-ray background. The integrated galaxy X-ray 
flux is given by 



4-kHq 



H (1 + • 



(l + z)(Q m (1 + Zf + OO 1 / 2 



dz, 



(8) 



where j x is the X-ray emissivity at z = 0, p is the adopted 
evolution and a x = T — 1, with V = 1.4. We integrate lumi- 
nosities above 10 41 ergs _1 assuming a pure luminosity evo- 
lution with p = 3. For the intensity of the XRB we adopt the 
normalisation of Vecchi et al. (1999). Integrating equation 8 
to 2 = 0.4 and z = 1 we find that the present sample con- 
tributes about 6 and 17 per cent respectively of the XRB. 
These fractions are somewhat uncertain depending on many 
assumptions regarding the evolution of the different X-ray 
populations and the absolute normalisation of the XRB in- 
tensity. We also note that the sample used in this calculation 
comprises both type-I and II Seyferts at moderate-z. There- 
fore, the XRB contributions estimated above include many 
of the nearby counterparts of the high-z type-I QSOs, which 
are responsible for a significant part of the XRB. 
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logL* gi 32* 4>+ jx 

(ergs -1 ) (Mpc~ 3 dcx -1 ) (erg s" 1 Mpc~ 3 ) 

43.37+° 0.72lg;^ 2.66±J° 52 1.03 x 10~ 5 1.9 x 10 38 



*The upper limit of the 32 index is unconstrained. See text for details. 
Table 4. Maximum likelihood parameters for the double-power law luminosity function of red X-ray sources with 2 < 0.4. 



5 DISCUSSION AND CONCLUSIONS 

We combine data from the XMM-Newton /2dF survey with 
public XMM-Newton fields overlapping with the SDSS to 
explore the nature of the hard (2-8 keV) X-ray selected pop- 
ulation hosted by red optically extended sources. The op- 
tical properties above suggest moderate redshifts (z < 1) 
and either dust reddening or optical continuum dominated 
by stellar emission. The significance of this X-ray popula- 
tion is threefold: Firstly, these sources comprise a signifi- 
cant fraction (about 75 per cent) of the moderate- z X-ray 
population to the flux limit of our survey, fx (2 — 8 keV) = 
2 x 10 -14 ergs -1 cm -2 . The X-ray luminosity function of 
the subsample with z < 0.4 is also found to be in fair agree- 
ment with previous estimates at low-z, further suggesting 
that our red colour selection picks the dominant X-ray pop- 
ulation at these redshifts. Secondly, using this luminosity 
function under the assumption of luminosity evolution of 
the form cx (1 + z) 3 we estimate that these systems could 
be responsible for a non-negligible component of the XRB 
to z = 1, about 17 per cent. Thirdly, we argue that our 
sample comprises a large fraction of the sources responsible 
for the spectral properties of the XRB. The stacked X-ray 
spectrum of the red X-ray sources is found to be consistent 
with a power-law with F ~ 1.4, i.e. similar to that of the 
XRB. 

Our sample comprises a mixed bag of objects including 
both absorbed and unabsorbed systems. Limited spectro- 
scopic information available for our sample suggests that 
the X-ray absorbed population is mostly associated with 
narrow-emission line galaxies (e.g. Seyfert-2 activity) while 
the unabsorbed sources show either broad emission-lines 
(e.g. Seyfert-ls) or absorption lines (e.g. XBONGs). We ar- 
gue that the red colours of the majority of the sources in the 
sample can be explained by the host galaxy stellar emission 
significantly contributing to the optical continuum. 

For the obscured X-ray population previous studies in- 
deed suggest that the optical colours of this class of sources 
are dominated by the host galaxy (Barger et al. 2002, 2003; 
Georgakakis et al. 2004a) with the AGN emission most likely 
blocked from view by the obscuring screen of gas and dust. 
In the case of the unobscured sources in our sample with 
Lx ~ 10 42 — 10 44 ergs -1 it is not clear why their colours arc 
redder that those of AGNs. For these systems we estimate 
a lower limit for the stellar contribution to the r-band op- 
tical continuum of 40 — 60 per cent. We show that an old 
(red) stellar population contributing up to 80 per cent in 
the r-band can indeed make these systems appear redder 
than typical AGNs. Higher fractions are however, required 
for young stellar populations. Severgnini et al. (2003) argue 
that the continuum emission of an Lx(2 — 8) = 10 43 ergs -1 
AGN can be swamped by a galaxy brighter than Mr = 
—22 mag. The median absolute magnitude of our sample is 



Mr — —22.32 mag, providing further support to the scenario 
of AGN dilution by stellar light. This luminosity is about 
1 mag brighter than the characteristic absolute magnitude 
of the luminosity function, M R — 21.2 (Blanton et al. 2003; 
Ho = 70kms -1 Mpc -1 ). Higher optical luminosities may be 
hard to understand in terms of stellar light alone. In these 
most extreme systems (in term of absolute magnitude) the 
AGN may dominate the optically luminosity and moderate 
dust obscuration (consistent with the estimated low Nh of 
the broad-line systems) may be responsible for the red opti- 
cal colours. For example in the case of Galactic dust-to-gas 
ratio a column density of Nh ~ 10 21 cm -2 corresponds to 
g — r reddening E(g — r) « 0.3 mag. 

Francis, Nelson & Cutri (2004) used the Two Micron All 
Sky Survey (2MASS) to explore the fraction of red AGNs 
at low-z likely to remain unidentified in UV/optical selected 
samples. They find that a small fraction of the J — K > 
1.2 mag selected population (about 1 per cent) is associated 
with broad emission-line AGNs. They also argue that host 
galaxy light contamination is responsible for the red colours 
of at least half of these systems. 

Cutri et al. (2001) also identified broad line AGNs (in- 
cluding types 1, 1.5 and 1.8) in the 2MASS catalogue albeit 
with extremely red colours (J — K > 2 mag), the vast ma- 
jority of which lie at z ~ 0.25. Wilkes et al. (2002) explored 
the X-ray properties of these broad line AGN and found ev- 
idence for obscuration in the range Nh ~ 10 21 — 10 23 cm -2 . 
They argue that dust reddening is most likely responsible 
for the red colours of these sources. However, the Wilkes et 
al. (2002) AGNs are, on average, strikingly underluminous 
at X-ray wavelengths relative to their near-IR flux. This 
is unlike the broad-line red sources studied here that have 
log fx /fopt typical of QSOs (see Figure 6). Moreover, the 
median Nh of the red X-ray sample studied here, including 
a substantial number of narrow-line systems with high ob- 
scurations, is « 10 21 cm -2 . This is lower than « 10 22 cm -2 
for the Wilkes et al. (2002) 2MASS AGNs, although the ma- 
jority of their systems have broad emission lines (types 1.8 or 
earlier) with only a small fraction showing narrow emission 
lines. Focusing on the spectroscopically identified broad-line 
sources in the present sample we find X-ray spectral prop- 
erties suggesting little photoelectric absorption above the 
Galactic value. We conclude that at least the spectroscopi- 
cally identified red broad-line AGNs in our sample are not 
the same population as those with very red optical/near-IR 
colours studied by Wilkes et al. (2002) but similar to the 
less extreme (in terms of J — K colour) red AGNs discussed 
by Francis et al. (2004). We cannot nevertheless, exclude 
the possibility that extreme objects like those studied by 
Wilkes et al. (2002) are present in our non-spectroscopically 
identified sample. 
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